Interplay of charge and heat transport in a nano-junction in the out-of-equihbrium 

cotunnehng regime 
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We study the charge transport and the heat transfer through a nano-junction composed of a small 
metallic grain weakly coupled to two metallic leads. We focus on the cotunneling regime out-of- 
equilibrium, where the bias voltage and the temperature gradient between the leads strongly drive 
electron and phonon degrees of freedom in the grain that in turn have a strong feedback on the 
transport through the grain. We derive and solve coupled kinetic equations for electron and phonon 
degrees of freedom in the grain. We obtain the heat fluxes between cotunneling electrons, bosonic 
electron-hole excitations in the grain, and phonons, and self-consistently find the current-voltage 
characteristics. We demonstrate that the transport in the nano-junction is very sensitive to the 
spectrum of the bosonic modes in the grain. 

PACS numbers: 72.15.Jf, 73.63.-b, 85.80.Fi 



I. INTRODUCTION 

Quantum composite materials are the subject of cur- 
rent intensive research activity. A typical quantum com- 
posite material (QCM) is a granular conductor, where 
the grains are so smah that the typical charging energy 
of one electron, E^., is the largest energy scale, in partic- 
ular it is larger than temperature and voltage. 

The equilibrium properties of granular conductors are 
well understood.lil However, much less is known about the 
transport properties of quantum composite conductors 
out-of-cquilibrium, e.g., at large bias or at large enough 
temperature gradients, when linear response theory is not 
valid. 

There are several transport channels for electron prop- 
agation though the granular QCM. At weak coupling 
between the grains and not very high temperatures cct 
tunneling is the main mechanism of electron transport .□ 
It provides a conduction channel at low applied biases 
and temperatures, where otherwise the Coulomb block- 
ade arising from the electron-electron repulsion would 
suppress the current flow. The essence of a cotunneling 
process is that an electron tunnels via.virtual states thus 
bypassing the large Coulomb barrier J3 Here we focus on 
the out-of-equilibrium cotunneling through a single grain, 
see Fig. |^), the building block of a quantum composite 
material, attached to two bulk leads. Voltage and/or 
temperature differences between the leads strongly drive 
electron and phonon degrees of freedom in the grain. The 
solution of the transport problem implies a selfconsistent 
calculation, which takes into account the mutual feed- 
back of highly excited degrees of freedom in the grain on 
the cotunneling electrons and vice-versa. 

Inelastic cotunneling of electrons from one lead to the 
other through a grain is accompanied by the creation of 
electron-hole pairs in the granule that results through its 



decay in heating of the grain. On the other hand, elec- 
trons in the grain exchange energy with phonons. Thus, 
the electron temperature in the granule depends on the 
balance between heating by electron-hole pairs and cool- 
ing by phonons. Recently heating effects in a single grain 
and a chain of grains due fjQ [inelastic cotunneling were 
discussed in several papers.oQ However, our considera- 
tion was in the linear response regime, i.e., valid for small 
gradients of voltage and temperature. Here we follow a 
general approach, which also holds in non-equilibrium 
situations and, furthermore, derive our results from mi- 
croscopic considerations. 

In the following sections, we first present our model 
and main results, followed by their derivations, and then 
discuss our results in detail. 



II. MODEL AND MAIN RESULTS 

We investigate out-of-equilibrium properties of a nano- 
junction shown in Fig. ^). The single grain in the junc- 
tion is characterized by two energy scales: (i) the mna: 
energy level spacing 5 and (ii) the charging energy EcU' 
Here we are interested in the case of metallic grains, 
where E^ 3> 5. In addition to theses two energy scales 
the system is characterized by the bare tunneling resis- 
tance i?L(R) between the left (right) lead and the grain. 
In this paper we concentrate on the Coulomb blockade 
regime, corresponding to a weak coupling between the 
grain and the leads, i?L(R) S> i?q with i?q = 7r?i/e^ being 
the quantum resistance. Fig. |l|. Since we focus on the 
cotunneling regime in this work, our considerations are 
valid for temperatures 5 < T < Ec and the leads are as- 
sumed to be heat sinks and stay at ambient temperatures 
Tl R, since they are much larger than the grain. 

As a main result, it turns out that the transport char- 
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FIG. 1. (Color online) Sketch of a nano-j unction, a), and its 
equivalent circuit scheme, b). 

acteristics of these junctions are very sensitive to the tem- 
perature of the electron degrees of freedom in the grain, 
Tg. The grain temperature, Tg(V, Ti^,Tji), is found exphc- 
itly as a function of voltage bias V and both lead tem- 
peratures from the expression for the heat fluxes in the 
nano-junction out-of-cquilibrium. Finally, this tempera- 
ture is used to obtain the current- voltage characteristics, 
I{V). 

The heat, Q, generated in the grain due to inelastic 
electron cotunneling can be derived as follows - repre- 
senting the first main result of our work: 

Q{T,,T^,,V) = a[T,\-T^], (1) 

Te'ff -Ti + ^TlieVr + ^ieV)\ (2) 

where = (T^ +T^) /2 and a material dependent con- 
stant a, which is derived below [see Eq. (|29|)1. The energy 
scale Teff has the physical meaning of an effective tem- 
perature of the cotunneling electrons. It is important to 
note that the r. h. s. in Eq. can be positive or negative 
indicating a heating or cooling of the grain, respectively. 
For equal temperatures Tl — Tr and zero voltage V ^ 
we reproduce the results of Ref. |^ from Eq. (|^). 

The heat dissipation in the granule goes through two 
stages. First, electron-hole pairs are excited by inelas- 
tic cotunneling processes and second the electron-hole 
pairs recombine and release their excess energy to the 
phonon bath. The heat flux between p hano ns and elec- 
trons (electron-hole pairs) has the formtH[ll3 

g(T„rph) = (T;-rp"j«, (3) 

where Tph is the temperature of the phonon bath and the 
parameters k and a depend on the particular model of 
electron-phonon interaction; typically a = 4, 5 or 6. 

In order to obtain the grain temperature, we have to 
solve the heat balance equation in the grain 

g(r<„rph) = Q(r„T„F), (4) 

where Q{T^,T^,V) is given in Eq. (|l|). 
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FIG. 2. (Color online) (top) Behavior of the grain tempera- 
ture Tg as functions of T^s- The blue curves are for a = 4, 
red for q = 5, and green for a = 6. The solid line are for 
Tph = 0.1, the dashed line for Tph = 0.5, and the dotted thick 
lines for Tph = 0.9. In all curves the dimensionless parame- 
ter is fixed to & = 0.5 (see text), (bottom) Behavior of Tg as 
functions of T^i for a = 4, Tph = 0.1, and different b. The 
solid curve is for 6 = 0.01, dashed for b = 1, dotted for b = 10, 
and dash-dotted for b — 100. The dependence of Tg on Tph 
for fixed Tes looks similar if b is replaced by b^^ (exactly the 
same for a = 4). 

It is useful to express all temperatures0 and voltages 
(more precisely eV) in terms of the charging energy Ec, 
which we will denote by a tilde, e.g., Tcff — Tcs/Ec or 
V = eV/E,. 

Introducing the dimensionless parameter b = 
aE^~°'/K, Eq. (^) can be written as 

Tg - Tp", = b (f^^ - f 4) . (5) 

For a = 4, Eq. (|) can be solved analytically: 

f^^{T^,, + bT,%)/{l + b). (6) 

The behavior of the grain temperature Tg as a function 
of Teff for different a- values and values of parameter b ~ 
0.01, 1, 10, 100 is presented in Fig. |. 



3 



Using Eq. (||) for the grain temperature Tg, we can 
finally calculate the current-voltage characteristics due 
to inelastic cotunneling 

(7) 

which is the second main results of this paper. 

When deriving Eqs. (|l|)-(||), we neglect possible se- 
quential electron tunneling. The latter is suppressed 
by the Coulomb blockade if temperature and voltage 
are below the characteristic single-electron charging en- 
ergy in the grain, eV,T < Ec with Ec ~ e^/ G^, where 
Ce = Cl + Cr + Cg is the total capacitance, see Fig. |l|b). 
Furthermore, the leads are assumed to be in equilibrium 
at similar temperatures and energy relaxation processes 
within the grain are fast such that we can use a local 



equilibrium description. 

In the following section we present the general out- 
of-equilibrium description of the heat and charge trans- 
port through the junction, which yields the derivation of 
Eqs. (|l|)-@ and eventually results in Eq. (0). 

III. COTUNNELING TRANSPORT 

We start with the most general expression for the 
current-voltage characteristics due to inelastic electron 
cotunneling. 

A. Inelastic cotunneling rate 

The inelastic cotunneling rate for an electron transfar 
from the left to the right lead has the following form:tll 



1 



1 



£2 - ei + AF+ (n) £4 - £3 + AF' (n) 



6{eV + ei~e2 + e3-e4), (8) 



where V is the bias voltage between the left and right 
leads. The distribution functions /^^'"^^(e) describe elec- 
trons within the left and right leads, respectively. 

If the temperature of the left (right) lead is T^(n) then 
/L(K)(e) = /F(e,T,(R)), with (e,T) = l/[exp(e/T) + 1] 
being the Fermi-Dirac distribution function. Function 
f^^\t) in Eq. (^) is the electron distribution in the 
grain. The functions l\F^^^-^{n) denote the change of 
the charging energy when an electron tunnels from the 
left (right) lead onto the island with excess charge ne, 
while AFj^j^j (n) denotes the change of the charging en- 
ergy when an electron tunnels out of the island with ex- 
cess charge ne to the left (right) lead 



e 



± n 



± n 



CgVg 



CaVa 



± 



[C„+Cgl2]V 



(9a) 



[C. + Cg/2]V 



(9b) 

where all the capacitances are defined in Fig. (|l)b) and 
Ce is the total capacitance introduced below Eq. (^). The 
number of excess electrons on the grain is determined by 
the condition that n results in the minimal electrostatic 
energy: Ai^+(n) > and /^F^{n) > 0,B 

The inelastic electron cotunneling rate from the 

right to the left lead can be found using Eq. (0) with 



proper substitutions: V — > —V in the 5-function and 
interchange between the left and the right leads, "L" ^ 
"R". 



B. Cotunneling current 

To write the total cotunneling current we also take 
into account elastic electron cotunneling rates, I^('=') and 
which can be written similarly to Eq. (H), i.e.. 



I{V) = e (l^('"'=') - t^(in°l) + r^(<=I) _ 



(10) 



Here we stress that the rates and the current in Eq. ( |lO| ) 
strongly depend on the electron distribution function 
/(9)(e) in the grain. This distribution can be approxi- 
mated by an equilibrium distribution with phonon bath 
temperature, Tph, only for small voltage and tempera- 
ture gradients across the junction. M In this paper we 
consider the current-voltage characteristics beyond the 
equilibrium approximation. We write the kinetic equa- 
tion governing the behavior of function f^^\e) and find 
the dependence of grain temperature Tg on the bias volt- 
age V, lead temperatures Tl, Tr, and the phonon tem- 
perature, Tph. 

In the following we concentrate on the inelastic com- 
ponent of the cotunneling current /('"°')(T/) which heats 
the grain in contrast to the elastic part. 



4 



IV. KINETIC EQUATION FOR /^^^ 

To calculate the current- voltage characteristics, I{V), 
in Eq. (|l^) we need to know the electron distribution 
function Z^^-* (e) in the grain. This distribution function 
satisfies the kinetic equation 



d{5t) 



where S is the mean level spacing in the grain. The 
left hand side of this equation describes the change of 
the electron distribution function with time t. The right 
hand side is the sum of collision integrals with ^('"'^'^(e) 
being the scattering integral due to cotunneling pro- 
cesses, while I*^''~P^^(e) and I'-'^^'^^e) being the electron- 
phonon and electron-electron scattering integrals, respec- 
tively. We write the scattering integral (e) explic- 
itly in Appendix |a[ 

We emphasize that only inelastic cotunneling con- 
tributes to the scattering integral l(™''')(e) in Eq. (11). 
It satisfies the particle conservation law 



X('"°')(e)de = 0. 



(12) 



This property does not contradict the current flow 
through the junction since the grain can be considered 
as a virtual haven for tunneling electrons during the co- 
tunneling processes. However, the scattering integral 
X(in<=i)(g) does not conserve energy, since each inelastic 
cotunneling process leaves an excited electron-hole pair 
behind in the grain. The heat dissipation rate into the 
grain is 



(13) 



Here the energy e in the electron grain distribution func- 
tion /^^^ (e) is counted from the local electrochemical po- 
tential. 



A. Local equilibrium approximation 

To solve the kinetic equation (^l]) we use a local equilib- 
rium approximation. This approximation is valid because 
for small grains the electron-electron (Coulomb) interac- 
tion is strong. Therefore the effective electron-electron 
scattering time, r^-e, corresponding to the scattering in- 
tegral I^'^~'^^(e) in Eq. ( pi] ) is shorter than the inelas- 

inoi and the electron- 



tic cotunneling scattering ti: 
phonon scattering time, To_ph,L. 



(14) 



Te-e < min(Tinol, To_ph). 

In addition, the scattering integral I^^""^' (e) in Eq. (Jll 
conserves the energy and the particle number. Therefore 
we can find the solution of the kinetic equation (|ll|) using 
the local equilibrium approximation 



This expression substituted into the scattering integral 
in Eq. (0) lets it vanish. Using the particle 
conservation law, Eq. (^2|), the correction to the effec- 
tive electro-chemical potential in the grain in the local 
equilibrium approximation can be neglected. Therefore 
in this approximation there is only one unknown param- 
eter, the grain temperature Tg. 



B. Heat balance equation 



The heat rate between electrons and phonons. 



(16) 



in the local equilibrium approximation, Eq. ( [l5| ) , is given 
by Eq. (|^). Therefore the problem of finding the grain 
temperature Tg in the local cqualibrium approximation, 
Eq. ([l5|), reduces to the solution of the heat balance equa- 
tion that follows from Eq. (ITT 



(17) 



q{Tg,Tpi,)^Q{Tg,T,,T^,V). 



Solving Eq. (0) we find the grain temperature 
TgiT^jTj,, V,Tph) as a function of lead temperatures Tl, 
Tr, the bias voltage V, and the phonon temperature Tph. 
Using the result for grain temperature Tg{Ti^, T^, V, Tph) 
we find the current- voltage characteristics in Eq. (px|). 
Below we proceed with this program focusing on the an- 
alytical solution of Eqs. ( p^ ) and (pT|). 



BOSONIC REPRESENTATION OF CHARGE 
AND HEAT RATES 



Here we consider the expressions for the inelastic co- 
tunneling rate r''"''') in Eq. (||) and for the heat dissipa- 
tion rate Q in Eq. (|l3|) in more details. 

For small grains the electrostatic energies, AFi^,~^ (n) 
and AFr,~ (n) in Eq. (^), are much larger than all other 
characteristic energy scales in the problem including tem- 
perature and voltage. Therefore we can neglect in the 
denominators of Eqs. (^), (Al), and (A3) the energy e 
and Ei, i = 1, ... ,4. This limit allows for an analytical 
solution of the heat balance equation ( p^ ) in order to cal- 



culate the grain temperature Tg. In addition, we show 
that it is more convenient to rewrite all transport char- 
acteristics in terms of electron-hole (dipole) excitations 
to find the current- voltage characteristics in Eq. (p^. 

First, we perform the substitution ei — > ei — eV/2 and 
€4 ^ ei + eV/2 in Eq. (||), which moves the voltage de- 
pendence into the distribution functions and in particu- 
lar implies that /L(R)(e) = /^^(e =F eF/2, rL(R)). After the 
(15) energy-integration, we obtain the following expression for 



5 



the inelastic cotunneling rate (see Appendix ^0 

^(incl) 



1 K K 

27r3 h ' 







dw ^ {nL^-> [N^ + !] + [! + nL^"']iV^} , (18a) 



where we used the notation 
1 _ 1 



1 



(18b) 



For the heat dissipation rate in the grain we obtain, using 
Eqs. (^)-(|a|): 



1 Rl 



7r3 i?Li?R hEl 

oj^duj {n^[N^ + l]-[l+n^]N^} . (19) 



Here we introduce the following form- factors, see Fi|g-^, 
that describe the nonequilibriuni electron-hole pairs 



,(LR) 



&/(^)(e+)fl-/(^)(e_) 



l + n^' 
1 + nr' 



— UJ 1 
(RL) 



(20a) 

(20b) 

(20c) 
(20d) 



where e± = e ± a;/2. General identities coupling the 
bosonic fo rm-f actors with and without tilde [see proof in 
Appendix B2] are given by: 



eV 
_eV_ 

U! 



(21) 
(22) 
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The functions and satisfy the basic property 
of the Bose-function Nb{uj,T) ^ l/[exp(w/r) - 1], 
Nb{-^,T) = -[1 + iVB(w,r)], e.g., = -[1 + n^]. 
Below we refer to the distribution functions and 
as to the electron-hole distribution functions in the leads 
and in the grain, respectively. 



A. Physical interpretation of equations (^) - (23) 



The effective electron-hole distribution function riu) in 
Eq. ( [23aD has a direct physical meaning. It describes the 
concentration of electron-hole pairs where the electron is 
localiaed-jn the left and the hole in the right lead or vice 
versalljo. 

The form- factor nJ_J'"' in Eq. (20a) is more specific 
than rii^. It describes the concentration of "polarized" 
electron-hole pairs with electron being sitting strictly in 
the left lead while the hole being occupying the right 
lead, see Fig. Therefore the annihilation of such an 
electron-hole pair (or the creation of the "RL" electron- 
hole pair) leads to charge transfer from the left to the 
right lead. 

Th e ex pression, n^i^^'>[Ni^ + 1], in the r. h. s. of, 
Eq. ( 18a ), is proportional to the probability of annihi- 
lation of an electron in the left lead and a hole in the 
right lead and the creation of an electron-hole pair lo- 
calized in the grain. The physical realization of such a 
process leads to a charge transfer from the left to the 
right lead. The product, [l-\-n'-^^'^\N^, is proportional to 
the probability of the creation of a hole in the left lead 
and an electron in the right lead and the annihilation of 
an electron-hole pair localized in the grain. The realiza- 
tion of this process also leads to a charge transfer from 
the left to the right lead. 

The expression, n^[N^ + l]j in the r. h. s. of Eq. ( |l9| ) 
is proportional to the probability for the annihilation of 
an electron-hole pair in the leads (electron in the left 
lead and a hole in the right lead and vice versa) and the 
creation of an electron-hole pair localized in the grain. 
This process increases the energy in the grain, but does 
not necessarily imply a charge transfer between the leads. 
Finally, the product [l-|-n^] A'^ corresponds to the reverse 
process that cools the grain. 



FIG. 3. The form-factor nj^^"^' (^L^"^') can be interpreted as 
the distribution function of electron-hole pairs where the elec- 
tron sits at the first lead while the hole on the second lead. 
The form-factor nJJ^'^' {fi'-^^'' ) corresponds to the opposite sit- 
uation. 

Finally we define the form-factors rii^ and A^^^ in 
Eqs. (|l|) 



1 



(RL)1 



(23a) 
(23b) 



B. Cotunneling current-voltage characteristics 



Using Eq. (18a) we can write the inelastic cotunneling 
current in the following form 



j(incl) 



I Rl e 



27r3 i?Li?R h J a E, 



duj — X 



nif-']A^^}. (24) 



Identities (^)-(|2^) allow us to simplify Eq. ( p4| ) and sep- 
arate the grain and the lead degrees of freedom. So we 
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get for the integral in Eq.(G4 



j(inol) 



1 



27r3 i?Li?R h 



duj 



2eV 



The integration over the energy uj in Eq. 
the final result for the inelastic cot unn e ling 
in Eq. (|o|), where we used Eqs. ( ^ ) 
gration of the n[ 



(25) 

(H) leads to 
eurrent given 
I for the inte- 



being the effeetive temperature, 
tion to the integral in Eq. (P7 



The second contribu- 
leads accordingly to 



uj^N^ duj = T^TT^/IS. Thus we obtain our main re- 
• suit, Eq. (hi), with parameter a given by the expression 



1 



15 i?Li?R hE^ ' 



(29) 



Here we 
Eq. (M, 



stress that the effective temperature Toff in 
follows the limit, lim„__j.o wn^^ 



eV 



eV 



2 COth 2^ 



(LR),(RL) 



parts of the integral in Eq. (p5t). only for smaU bias vohages, eV < T„ 



C. Evaluation of the heat dissipation rate Q in the 

grain 

Formally the heat dissipation rate into the grain Q in 
Eqs. (|l3| ) and (^9|) can be obtained using the effective 
kinetic equation for electron-hole pairs 



d 



d{h-^E^t) 



r(incl)(-^) 



(26a) 



where the relaxation rate for electron-hole pairs in the 
grain is mediated by their interaction with cotunneling 
electrons 



1 



1 R, 



2 



(26b) 



Equation (26b) means that the relaxation rate of 
electron-hole pairs in the grain in the cotunneling chan- 
nel is proportional to the second power of temperature 
at small temperatures, 1/t('"°')(w) ~ T^. 

For very low temperatures, where electron-electron and 
electron-phonon interactions are frozen out and scatter- 
ing due to cotunneling is the leadin g sc attering mecha- 
nism we obtain ~ n^^ using Eq. (26a). 

However, for a small metallic grain the electron- 
electron interaction is the main scattering mechanism. It 
drives the electron- hole distribution function, N^: , to the 
local-equilibrium form, Ni^ « Nb{ui, Tg). [The last state- 
ment follows from the fact that in the local equilibrium 
approximation the distribution function /^^^ « fp{t,Tg) 
and therefore according to Eq. ( ^3b| ), N^^ w NB{uj,Tg).] 
In this case the effective grain temperature Tg is deter- 
mined by the heat-balance equation (^7|). 

To evaluate the heat dissipation rate in the grain, Q in 
Eq. (|l9|), we need to know the integral 



Q 



1 Rl 



1 



7r3 i?Li?R hEl 



{nu - N^) duj. 



(27) 



Using Eq. (|23al) and Eqs. (p5|)-(p8|) (see Appendix H) for 
ji m 
where 



the first term in this integral we obtain uj^n,,, duo 



r„1 



rj-rl 



(28) 



VI. DISCUSSION 



Quasi-equilibrium limit of the current voltage 
charateristics 



For equal lead and grain temperatures, Tl = Tr = 
Tg, we reproduce the known result for the cotunneling 
current, Refs. [|l| and O. However, in the general case, 
the inelastic current in Eqs. ( p4| ) and (|^) strongly 

depends on both the grain and the lead temperatures. 



B. Heat dissipation rate and total power 

The total power of the junction, given by the product 
of current and bias voltage, / V and the heat dissipated 
into the grain, Q{Tg,Ti,T2, V), are two different quanti- 
ties. Indeed, the heat Q{Tg,Ti^,T„,V) is finite even for 
zero bias voltage, V = 0, ii the lead temperatures are 
different from the grain temperature, Tl^r 7^ Tg, while 
the total power is zero, I V ~ 0, in this case. 

For finite voltage, V ^ 0, and equal lead temperatures, 
Tl = Tr, the heat dissipation rate Q can be identified 
with the Joule heat released in the grain. There is an- 
other part of the Joule heat, Qioads, corresponding to the 
energy equilibration of nonequilibrium electrons in the 
bulk of the lead after they co-tunnel from the other lead. 
Therefore the Joule heat is the sum of both contributions, 

IV ^ Qlcads + Q- 

It is interesting to compare the heat dissipation rate 
in the grain Q with the Joule heat, IV. Using Eq. (|l|) 
for the heat dissipation rate in the grain Q and an ex- 
plicit expression for inelastic part of cotunneling current 
Eq. (0) we obtain 



47r4 



rp4 



a 



5 {(eV^)2 + 27r2(T2 -HT2)} (eF)2 



(30) 



First, we consider the case of equal temperatures, T = 



Ti^ = Tji = T„. In this case we obtain 



Q 



/(incl) Y 



(31) 



Equation (^Tj) has transparent physical meaning - half of 
the energy in the cotunneling process is spent to generate 
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electron- hole pairs in the grain and the other half to inject 
the nonequilibrium electrons into the leads that finally 
equilibrate in the bulk of the leads. Remarkably, the 
ratio F) is universal and does not depend on 

temperature T and voltage V . 



C. Transport in semiconducting nanojunction and 
spectrum of bosonic modes in the grain 

Calculating the current and the heat we implicitly as- 
sume that the grain is metallic meaning that electron- 
hole pairs with an arbitrary small energy lo can be cre- 
ated during the inelastic cotunneling process. However, 
if the grain is made of a semiconducting material the 
spectrum of electron- hole pairs would have a gap A. In 
this case the lower limit of integration over the energy 
uj in Eq. ( p^ ) would be A and a smooth weight func- 
tion p{u)) would appear under the integral renormalizing 
the interaction vertex of electron-hole pairs with the co- 
tunneling electrons. Therefore in this case the current 
would be exponentially suppressed for voltages V and 
temperatures T ~ Ti^ = T„ = Tg smaller than the gap A, 
^ y exp(-A/r) . This shows that the current- 
voltage characteristics is very sensitive to the spectrum 
of the grain. 



VII. CONCLUSIONS 

We studied the charge transport and the heat trans- 
fer through a small metallic grain weakly coupled to two 
metallic leads. We focused on the cotunneling regime 
out-of-equilibrium, when the bias voltage and the tem- 
perature gradient between the leads strongly drive elec- 
tron and phonon degrees of freedom in the granule that 
in turn have a strong feedback on transport through 
the grain. We derived and solved the coupled kinetic 
equations for electron and phonon degrees of freedom 
in the grain, found the heat fluxes between cotunneling 
electrons, bosonic electron-hole excitations in the grain, 
and phonons. and selfconsistently obtained the current- 
voltage characteristics. We demonstrated that the trans- 
port in the nanojunction is very sensitive to the spectrum 
of the bosonic modes in the grain. 
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Appendix A: Scattering integral due to cotunneling 

Here we present the scattering integral l('"°'^(e) in 
Eq. (O) induced by cotunneling processes 



27r3 i?Li?R h 



1 1 



e-ei + AF+(n) e4-e3 + AF^{n 



S{eV + ei — e + €3 - 64)- 



1 R, 1 

27r3 R^Rj, h 



^ yiL, C4 - C3 -r L^± \ 

■ dei&2de4./^^'(ei)[l - /(^'(e2)]/'^'(e)[l - /'^'(e4)]x 



1 1 



. £2 - ei + AFl{n) £4 - e + AFj^ (n) 

where the last term in the brackets results from the contribution of the electron cotunneling rate from the 

right to the left leads. Elastic cotunneling does not contribute to the scattering integral. 

Using Eq. ( Al) we rewrite the heat flux dissipated into the grain in terms of the heat rates similar to Eq. (^: 

Q = fQ + fQ, (A2) 
where the heat rate from the left to the right lead is 

- ^j§^Ji I de.de^de^de^ie^ - e^) f'^'^\e,)[l - /(^)(e2)]/^^' (e3)[l - f'^"\e,)]x 

<5(eU + ei -e2 + e3-e4). (A3) 



t2-ti+ AF/ (n) £4 - £3 + AFfl {n) 

I 

The plus sign in Eq. ( [A2| ) implies that the particular right or from the right to the left, is not important for 
direction of the cotunneling process, from the left to the 
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heating or freezing of the grain. 



Appendix B: Bosonic representation of charge and 
heat rates 



1. Inelastic scattering rate 

Integrating the expression for the inelastic cotunnehng 
rate r (in^O in Eq. (|) over the E = {e2 + £3)72, E' = 
(ei + £4)72 and introducing new variables uj — £2 ~ ^3 
and cli' = £1 — £4, we obtain 

Jei,...,e4, 

J du;du;'dEdE'fl,';^[l - f^JP][l - fi'%'!d{Lo - J) - 

do.w^n[,^"'[l + 7V„], (Bl) 



where E± = E ±u/2 and E'+ E' ± u' /2. The form 
factors n'J'^' and in Eq. (Bl) can be interpreted as 
the effective distribu tion s of electron-hole pairs, see dis- 
cussions below Eqs. ( ^Oa| ) and (23a). This interpretation 
is possible for positive frequencies w > 0. Therefore it 
is co nve nient to transform the integral in the last line of 
Eq. (|l|) and find 



/ duji^^ [1 + N^] + n'Z' [1 + N^^] ] = 

Jo 

']N.}- (B2) 



/ 



[1 + n[^^>l 



This result immediately leads to Eq. ( |18a| ) for the inelas- 
tic scattering rate given in the text. 



2. Bosonic form- factors 



To pr ove identities, Eqs. (^l|)-(p2|), we rewrite 
Eq. (p^), 



1 + h; 



def^\e^) h -/(^•)(e+)' 



t(RL) 



UJ 



de 



/(«)(£+)- /(^)(£_)1 . (B3) 



Taking the distribution functions f^^\e,Ti) = fp{e ~ 
eV/2,Ti) and /(«'(£, T2) = /f(£ + 6^/2, T2) and using, 
JZoiM^'Ti) - Me,T2)]de = 0, wc find for the last 
integral: 



de 



de 



d 1 r 



uj + eV 



de2 



■ X 

eV 

UJ 



(B4) 



This proves the validity of Eqs. (pl|)-(p2|). 

Using the distribution functions f^^{e,Ti^) = fr je — 
eV/2,T) and f^^^K^^T^) = /f(£ + eV/2,T) in Eq. (P), 
with T being the temperature oLboth leads, we find the 
explicit form of the form-factoral3 



u — eV 



NBiuJ~eV,T), 



eV uj + eV 



UJ UJ 

LU + eV 



NB{LJ + eV,T), 



UJ 

eV 



NB{uj + eV,T), 



UJ — eV 



NB{Lo-eV,T). 



(B5) 
(B6) 
(B7) 
(B8) 



Equations (B5)-(BS) are exact for equal lead temper- 
atures only, T = Tl = Tr. However, if the lead tempera- 
tures are different ^ then Eqs. (^-(^ give still 
a good approximation of the form-factors if we use the 

substitution, T — ^ = y ^ ^ ^ . This is a good ap- 
proximation for Tl^r > eV and for Tl.r < eV when Tr ~ 
Tl. Furthermore, it reproduces lim„_j.o wn'J''^''"^^* = 

2~ COth 2^2^ ■ 

The integrals of form uj°'ni^d uj, a = 2, 3 , . . . c an b e 
calculated analytically using Eqs. (23a), (pa) and (B7), 
and expressed through the polylogarithmsEI So, for in- 
stance. 



uj^riu, duj = 



{ [Li3(e"") - Li3(e")] v + 3 [Li4(e-") + Li4(e^)] } 



-— ) 




V ^ 


27ri/ 


27ri/ 


(7rT)4 


\ 27r2 


5 


15 


8^ 



where we use the notation v = eV/T. 
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